On the effective size of certain "Schrodinger cat" like states 
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Several experiments and experimental proposals for the production of macroscopic superpositions 
naturally lead to states of the general form \0i)® N + \4>2)® N , where the number of subsystems N is 
very large, but the states of the individual subsystems have large overlap, |(<^i|<^2}| 2 = 1 — e 2 . We 
propose two different methods for assigning an effective particle number to such states, using ideal 
Greenberger-Horne-Zeilinger (GHZ)- states of the form 10)®" + |l)® n as a standard of comparison. 
The two methods are based on decoherence and on a distillation protocol respectively. Both lead 
to an effective size n of the order of Ne 2 . 

PACS numbers: 03.65.Ud, 03.67.-a, 03.65.Yz 



It was pointed out already in 1935 by Schrodinger 
that quantum mechanics predicts the existence of super- 
positions of macroscopically distinct states. The obser- 
vation of the related quantum interference effects is very 
difficult because of environment-induced decoherence. 
Nevertheless several methods for producing and veri- 
fying macroscopic superpositions have been proposed, 
in systems ranging from superconductors over Bose- 
Einstein condensates (BECs) || || and opto-mechanical 
systems || to small cantilevers coupled to supercon- 
ducting islands |J. Recently there have even been the 
first experimental demonstrations of the superposition of 
distinct macroscopic current states in superconducting 
quantum interference devices (SQUIDs) 0. 

The states produced in such proposed experiments can 
often be described to a good approximation by 



(i) 



with K = 2 + <0i |0 2 > + (</>2 1 01 ) ■ Here \ip) is a state 
of N two-level systems (qubits). The individual qubits 
could be seen as simple models for many different phys- 
ical systems, including the atoms in a BEC inside a 
double- well potential |3J , atoms in two internal states or 
the Cooper pairs in a SQUID (which can flow in clockwise 
or anti-clockwise direction). 

The essential point for our present discussion is that 
the states |0i) and |02) are not necessarily orthogonal. 
In fact, we will study the case where | (0i \4> 2 ) | 2 = 1 — e 2 is 
very close to one (e is small). Note that in spite of this the 
overlap between the two terms in ([l]) can be very small 
for large N, since it is given by |(0i|02)| 2Ar = (1 — £ 2 ) N , 
which is well approximated by e~ Ne for small e. 

We investigate how states of the form (Q) compare to 
ideal GHZ states of the form 



\GHZ n ) = ^=(\0)® n + \l)® n ). 



(2) 



States of the form (g) can be produced to good approx- 
imation in quantum optical systems, including atoms in 



cavity quantum electrodynamics Q , trapped ions II and 
photons from parametric down-conversion JTo| | . So far, 
particle numbers n up to 4 have been achieved. On the 
other hand, states of the form (|]) can involve macro- 
scopic numbers of particles, in combination with small 
values of e. An important question is whether there is 
a well-defined way to compare these two - very differ- 
ent - cases. That is, can one give a meaningful answer 
to the question whether e.g. a state of the form ([[]) with 
N = 10 6 particles, but e = 10~ 3 is more or less entangled 
than an ideal GHZ state with n = 10? 

A first simple way of assessing the "size" of states of the 
form (]l|) is to look at the overlap between the two terms. 
However, as pointed out above, this will be close to zero 
in most interesting cases and thus does not lead to a very 
sensitive criterion. A very intuitive way of comparing (|l|) 
and (||) would be to assign to (|l|) an "effective particle 
number" n, which could be interpreted as saying that the 
state ([!]) is (in a certain well-defined sense) equivalent to 
an ideal GHZ state of n qubits. This requires well-defined 
and physically meaningful methods of determining such 
an effective particle number. Here we will propose two 
such methods and show that they lead to essentially the 
same result, namely that the effective n of a state of the 
form ([!]) is of the order of Ne 2 . 

Our two methods are very different. The first one is 
based on the rate of decoherence. An important potential 
application of states of the form ([!]) is for the observa- 
tion of (or the search for) weak decoherence processes, 
ultimately including unconventional ones as predicted by 
spontaneous wave-function collapse models Jll]]. If the 
state (Q) is as sensitive to decoherence as an ideal n-qubit 
GHZ state, it is natural to say that its effective size is n. 

The second method of assigning an effective size to our 
states is in the spirit of quantum information, viewing 
multi-party entanglement as a convertible resource. We 
ask how much ideal GHZ entanglement can be distilled 
from the states (Q) by local operations (acting separately 
on each qubit) and classical communication. 

It is remarkable that these two very different ap- 
proaches lead to the same result. We believe that this 



2 



suggests that Ne 2 is indeed a good physical quantifica- 
tion of the "size" of macroscopic superpositions of the 
present type. We conclude by giving another simple ar- 
gument in favor of the proposed scaling of the effective 
size, based on particle loss. 

Let us now follow our first approach and study the 
effect of local decoherence on the state (|l|). We will con- 
sider the case where each of the particles undergoes an 
independent decoherence process, i.e. each particle is 
coupled to an independent bath. The effect of decoher- 
ence is quantified as the rate of decay of the off-diagonal 
elements in the natural basis. Note that we are interested 
in properties of states and not of physical set-ups which 
generate them. In this sense, although for different phys- 
ical systems the decoherence process may be completely 
different, we can study the behavior of the states describ- 
ing those systems under a certain decoherence process in 
order to compare the properties of these states. 

We consider phase decoherence in the natural basis, 
that is diagonal elements in the { |0) , 1 1) } basis remain un- 
changed, while off diagonal elements |0)(1|, |1)(0| decay 
with a rate 7. The decoherence process of an individual 
system is described by — > e _7 *|?)(j| for i ^ j and 
— > which corresponds — in a quantum infor- 

mation language — to a dephasing channel. The action 
of this channel may be described by the completely pos- 
itive map £ defined through £{p) = Pop + (1 — Po)o z pa z , 
where po = (1 + e _7 *)/2 and a z is a Pauli matrix. 

It is straightforward to establish the effect of this deco- 
herence process on an ideal n-particle GHZ-state. The 
density matrix for the state Eq. (§) is l/2(|0)(0|® n + 
|l)(l|®n + |o)(l|®« + |l)(0|® n ). Since £® n (a® n ) = 
(£ (<r))® n , one can study the decay of the off-diagonal el- 
ements |0)(l|® n and |1)(0|®" by considering the action of 
the decoherence channel £ on the single-qubit operators 
|0)(1| and |1)(0|. To be specific, we consider the trace 
norm, = trV A* A, of a t — £(aa), where a = |0)(1|. 

Note that P 0JV ||i = ||A||f . Since a t = e~ 7t a , we have 
that ||af"||i = ||a t ||? = e- 7 ™*||a ||i. The off-diagonal 
element of the GHZ state decays with a rate jn. 

We want to compare this to the decay rate of the off- 
diagonal terms for the iV-particle states \ip) of the form 
(|l|) . Without loss of generality we denote 

I'M = |o>, 

|^a) = cos(e)|0)+sin(e)|f), (3) 

and use the shorthand notation c e = cos(e) and s e = 
sin(e). We have that K = 2(1 + c^) and for small 
£) I (0i 102) | 2 = c? ~ I — e 2 . The density matrix for 
the state Eq. (ffl) is l/KQ^)^]^ + l&X&l®* + 
\<Pi}(fa\® N + \<f>2)(<t>i\® N )- We are interested in the de- 
cay rate of the off-diagonal elements. As before, the 
problem can be reduced to studying single-qubit oper- 
ators, namely b = \(/>i){<fa\ = c e |Q> <0j + s e |0)(l|. We 
have that 60 changes due to the above decoherence pro- 
cess to b t = £(b ) = Vd\0)(xt\ with \ Xt ) = I/v^qIO) + 
s e e~ 7 *|I)) and d — c 2 + s 2 e -27 *, such that \\t) is prop- 
erly normalized. It follows that H&f^Hi = d N l 2 . For 



e < 1, t < 7" 1 , Ne 2 ^t < 1, we have d w 1 - 2e 2 ^t and 
thus d N / 2 » e^ N(2t . This implies that the rate with 
which the coherences of the state decay is given by 
jNe 2 . That is, the decoherence rate of a state of the 
form Eq. ([j]) with |(</>i|</>2)| 2 = 1 — e 2 is the same as that 
of an ideal n-party GHZ state with n — Ne 2 and thus 
one may associate an effective particle number n — Ne 2 
to the state \ip). 

The observed decoherence rate is not restricted to this 
specific decoherence model. Consider for example the 
basis independent decoherence model of a partially de- 
polarizing channel. In this case, the individual deco- 
herence process for each qubit is described by —> 
MN)01 + (1 ~ f^^idh^ where p = e~ 7t . Equivalently, 
the completely positive map £ describing this process is 
given by £ (p) = YA=oPi a iP a i witn Po = (3m + l)/4 and 
Pi = P2 = P3 = (f — A t )/4, where er — 1, and the Oi are 
Pauli matrices. We find that ao [bo] changes due to this 
decoherence process to a t — £(o,q) = pao [bt — £(bo) = 
c e (l+/x)/2|0)(0| + (l-/x)/2|l)(l|+s e /i|0)(l|]. One obtains 
that 1 1 a* 1 1 1_ = p = e~ 7t and ||6t||i = \/c 2 + p 2 s 2 = V", 
which is exactly the same as in the case of the dephasing 
channel. One thus recovers exactly the same decoherence 
rates -777. and 77Ve 2 respectively- as in the case of the 
dephasing channel. 

Let us now turn to our second approach, which is more 
in the spirit of quantum information. We again consider 
states of the form ([!]) with 1 0i ,2) defined in Eq. (|^). 
We are interested in the distillation of ideal n-particle 
GHZ states (||) from these states under the condition 
that only local operations and classical communication 
are allowed. The restriction to local operations is essen- 
tial if one wants to quantify the entanglement contained 
in a given state because non-local operations could create 
additional entanglement. We are only interested in the 
number of particles which form a GHZ state after the dis- 
tillation process, i.e. the effective size of the GHZ-state, 
and not which of the particles are entangled. 

We show that the average number of the particles 
which is in an ideal GHZ-state after the distillation pro- 
cess scales essentially like n = Ne 2 . We (i) provide an 
explicit protocol to produce — with unit probability — 
n-party GHZ states from a single copy of where the 
average value of n is Ne 2 /2 and (ii) show that even in 
the asymptotic limit, i.e. considering several identical 
copies of the state \ip) , this average value is bounded 
from above by n w Ne 2 (- log 2 (e)/2) @. 

Let us start with (i), a practical protocol which trans- 
forms a single copy of \ip) deterministically into n-party 
GHZ states by means of local filtering measurements. 
The protocol we propose is a generalization to multipar- 
tite systems of the distillation protocol of Ref. jl3| for the 
optimal distillation of tripartite GHZ states from a single 
copy of an arbitrary pure state of three qubits. Consider 
the local filtering measurement described by the oper- 
ator A = fc(|O)(0i| + |1)(0 2 |), where {|&},|^}} is the 
biorthonormal basis to {|</>i), |<^2)}, i-e. (<f>j\(/>i) = 5ji- 
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The constant k is chosen such that the other opera- 
tor A of the local, two-outcome generalized measure- 
ment {A, A} —which fulfills A^A + A^ A = 1— has 
rank one. This implies that in case one obtains the 
outcome corresponding to A, then A ® ~i® N ~ 1 \ r ^>) = 
fe/V^dO)^!)®^- 1 + ll)!^)®^" 1 ), while for the other 
outcome A® l®*" 1 |V>) oc \ X ) ® (l^i)^" 1 + l-fe)®^ 1 ), 
i.e. the measured particle factors out. The distillation 
protocol works as follows: Each of the parties performs 
locally the two-outcome generalized measurement {A, A} 
and all those n parties which obtained a positive out- 
come, i.e. the outcome corresponding to A, finally share 
an ideal n-party GHZ-state of the form (||), while the 
remaining parties are in a product state. We shall be 
interested in the expectation value of n. 

Using the notation of Eqs. (Q, |3|) , we find that 



.4 



I 



(4) 



In the j th measurement, the probability to obtain the 
outcome corresponding to A is given by p = (I — 
c e )/(l + c^ - ^ 1 ) provided that none of the previous 
measurements was successful. In case one of the pre- 
vious measurements was already successful, the proba- 
bility to obtain an outcome corresponding to A is given 
by p = (1 — c e ) for the remaining parties. This different 
behavior after the first successful measurement can be 
easily understood by noting that once one of the mea- 
surements was successful, then the (normalized) state 
after the measurement is given by l/\/2(\0}\<j)i}® N + 
1 1 } I </>2 } ® ^ ), while otherwise the normalization constant 
is K = 2(1 + c^-J +1 ). One finds that the probability 
q n to obtain n (where n > 1) successful measurements 
— and thus n-party GHZ states — is given by 



q n = (i 



Ce)"cf- 



N 



1 



1 



(5) 



while qo — 2cf/(l + ). Note that this probability 
distribution is — up to the factor 1/(1 + c£) and corre- 
spondingly the value of go — very similar to a binomial 
distribution. The expectation value < n >= Y^,jLo 1j3 ^ s 
given by < n >= (I — c e )N/ (1 + ) which simplifies for 
Kl and Ne 2 > I to 



< n >« Ne 2 /2. 



(6) 



This provides the desired lower bound for the distillation 
rates of n-party GHZ-states. 

Regarding (ii) , the announced upper bound for the dis- 
tillation rate, we use the fact that the von-Neumann en- 
tropy of the reduced density operator with respect to 
system I p\ |l4|], S\(pi) = — tr(pi log 2 pi), is an entan- 
glement monotone, i.e. not increasing under local opera- 
tions and classical communication |l^, [if], [fTj . 

We consider the distillation process in the asymptotic 
limit, i.e. the transformation of M — > oo identical copies 



of the state \ip) to n-particle GHZ states Such a dis- 
tillation protocol consists of an arbitrary sequence of local 
operations (including measurements), possibly assisted 
by classical communication, mathematically described by 
a multi-local superoperator |l6). The protocol produces 
a certain number, say M n copies, of n-party GHZ-states, 
which can — as we are only interested in the number of 
parties which constitute a GHZ state — without loss of 
generality be considered to be symmetrically distributed 
among the N parties JlS}]. Such a symmetric configu- 
ration is denoted by \GHZ n )® M ™ |§. The distillation 
protocol is described by the following transformation 



N 



\<g>M 



\GHZ 7 



(7) 



n=2 



where p n > denotes the average number of n-party 
GHZ states which are produced per copy from 

Given the monotonicity of the entropy under local op- 
erations, we obtain 



N 



MSM) > > >i(|Gi/Z„}^" M ) 



(8) 



n=2 



where Si(\GH Z n )®' PnM ) denotes the entropy of the re- 
duced density operator with respect to system I of 
p n M copies of (symmetrically distributed) n-particle 
GHZ states [fijj. Since the probability that the first 
particle belongs to the n entangled particles is given 

by p = ( m \]l\ m ) = n /^' we nave that 



5i(|GffZ„)® p " M ) = PnMn/N and thus E^= 2 P» n < 
NSi(\ip)). It is straightforward to calculate SiQip)) using 
that for \ip) ([]]), the reduced density operator with respect 
to system! is given by pi = [(I+c 2 +2cf )|0)(0|+s e c e (l + 



C f- 2 )(|0)(I| + |f)(0|)+^|l)(l|]l/( 2 + 2 C , 



For e < I 



and Ne 2 > 1, one obtains that Si(\ip}) « -e 2 log 2 (e)/2 
which implies 



N 1 
Y,Pnn<-Ne 2 -\og 2 (e) 



(9) 



71=2 



as announced p0| . 

We would finally like to mention another simple ar- 
gument that suggests the same scaling for the effective 
size of the states Eq. ([[]). Let us compare the effects 
of particle loss on the state ([!]) and on an ideal GHZ- 
state. Suppose that every qubit is lost with a probability 
A. Consider an n-qubit GHZ state. As soon as a sin- 
gle qubit is lost, the state becomes completely diagonal. 
There is only an off-diagonal element in the case of no 
losses, which has a probability of (I — A)™. The expecta- 
tion value of the off-diagonal element in the case of losses 
is therefore |(1 — A)", equal to \e~ Xn for small A. 

On the other hand, for the state ([!]), tracing out 
particles reduces the size of the off-diagonal terms but 
does not completely remove them. Tracing out k par- 
ticles multiplies the off-diagonal terms by a factor of 
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(0i|02> fc = (1 - e 2 /2) fc , equal to e~ ke ' V 2 for small e. The 
typical number of particles lost will be NX, therefore the 
typical off-diagonal term will go like e~ XN<i I 2 . For large 
NX the probability distribution will be strongly peaked 
around the typical value. Therefore the expectation value 
of the off-diagonal term will be of order e~ XN<i I 2 . One 
sees that if the ideal GHZ state has n = Ne 2 /2, then the 
expectation values of the off-diagonal terms will have the 
same size for the two states. This is one more confirma- 
tion for our proposal that the "effective size" of the state 
(Q) scales like Ne 2 . Note that there are several other ar- 
guments which confirm this effective size n w Ne 2 , e.g. 
an argument related to the statistical distinguishability 
of states |0i,2} as pointed out in Ref. 

To summarize, we provided two different methods to 
assign an effective particle number to GHZ-like states of 
the form |^) cx \<f>i)® N + \fa)® N with |(</>i|0 2 )| 2 = 1-e 2 . 
The first method is based on the rate of decoherence, 
and we found that \ip) behaves like an ideal n-party GHZ 
state with n w Ne 2 . In the second method, which is more 
in the spirit of quantum information, we provided lower 
and upper bounds for the distillation rates p n of ideal n- 
party GHZ-states using only local operations and classi- 
cal communication. Again, we found that ^2p n n rs Ne 2 , 
i.e. the average number of particles which form an ideal 
n-party GHZ-state essentially scales like n « Ne 2 . This 
illustrates that not only the number of particles but also 
the properties of the states appearing in the microscopic 
description of the system determine the effective size of 



the corresponding Schrodinger cat like state. 

Some open questions remain. On the one hand, we 
have considered a particular class of Schrodinger cat like 
states (0). One can have physical systems where the 
macroscopic superposition cannot be described by states 
of this form. In particular, the states appearing in the 
superposition may not be a tensor product of identical 
microscopic states, \<f>i,2) , either because they are en- 
tangled themselves — e.g. the position states of the atoms 
in an oscillating micro-mechanical cantilever or mirror — 
or the the corresponding system cannot be decomposed 
in a natural way into subsystems — e.g. a superposition 
of two coherent states, |a) + | — a). It would be interest- 
ing to study by similar means the degree of " catness" for 
those systems and maybe compare them with the ones 
treated here. 

On the other hand, there are experiments where 
macroscopic superpositions have been created but there 
is no microscopic description of the states produced . 
It would be interesting to find such a microscopic de- 
scription and — in case the states can be written as (|l|) 
(which would be natural) — to assess an effective size of 
the states for these experiments. 
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